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Abstract 

We study the evolution of non-linear spherically symmetric inhomogeneities in string cosmology. Fried- 
mann solutions of different spatial curvature are matched to produce solutions which describe the evo- 
lution of non-linear density and curvature inhomogeneities. The evolution of bound and unbound in- 
homogeneities are studied. The problem of primordial black hole formation is discussed in the string 
cosmological context and the pattern of evolution is determined in the pre- and post-big-bang phases of 
evolution. 

PACS numbers: 98.80Hw, 04.50.+h, 11.25.Mj 

1 Introduction 

Considerable attention has been devoted to exploring the range of behaviours displayed by the equations of 
string cosmology which are supplied by the variation of the low-energy effective action of the bosonic sector 
of string theory Investigations have been made into the evolution of isotropic cosmologies [0], simple 
homogeneous anisotropic cosmologies of Bianchi type, 0, and Kantowski-Sachs type §, by various authors, 
and the present authors have provided a systematic classification of spatially homogeneous string cosmologies 
in terms of their relative generality when considered as constrained systems of non-linear ordinary differential 
equations ||. 

As in the case of general-relativistic cosmologies, the introduction of inhomogeneities into the string 
cosmological equations produces a considerable increase in mathematical difficulty: non-linear partial dif- 
ferential equations must now be solved. In practice, this means that we must proceed either by means of 
approximations which render the non-linearities tractable, or we must introduce particular symmetries into 
the metric of space-time in order to reduce the number of degrees of freedom which the inhomogeneities can 
exploit. Accordingly, inhomogeneous string cosmologies have been investigated in the approximation of small 
perturbations of the isotropic Friedmann-Robertson- Walker (FRW) models, in the 'velocity-dominated' ap- 
proximation, and by studies of exact inhomogeneous solutions with cylindrical symmetry. Barrow and Kunze 
H found a wide class of exact cylindrically symmetric flat and open inhomogeneous universes. Closed cylin- 
drical solutions were then found by Feinstein et al. pj. These solutions provide exact 'velocity-dominated' 
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solutions of general relativity and are expected to form a leading-order approximation to part of the general 
solution of string cosmology in the neighbourhood of the singularity (if we ignore higher-order string correc- 
tions to the action). Veneziano |8| and Buonanno et al || have studied the behaviour of string cosmologies 
at early times in the velocity-dominated approximation, (that is, neglecting spatial derivatives, relativistic 
motions, and 3-curvature inhomogeneities with respect to time-derivatives in the field equations pQ| ). The 
asymptotic forms obtained at early times also approximate the behaviour displayed by the exact solutions of 
Barrow and Kunze on scales larger than the horizon where the inhomogeneities evolve slowly: the inhomo- 
geneities are just homogeneously propagated. However, the exact solutions also provide information about 
the evolution of dilaton, axion, and gravitational wave inhomogeneities after they enter the horizon, where 
they attenuate by non-linear oscillations because of the pressure forces exerted by the dilaton and axion 
fields. These solutions do not contain trapped surfaces and so they cannot be used to follow the collapse 
of inhomogeneities to black holes, although it would be possible to study this problem by using the closed 
(S 3 ) solutions studied by Feinstein et al in [Qj. 

The exact solutions given in refs. Q and Q possess cylindrical symmetry and all physical quantities 
depend on at most one space coordinate and the time. The case of cylindrical symmetry is natural because 
of the mathematical simplicity of the field equations whenever there exists a direction in which the pressure 
equals the energy density. However, it is also important to consider the case where the inhomogeneities 
possess spherical symmetry. Not only does this seem more natural, in that there need exist no preferred 
direction in which the inhomogeneity dominates, but it allows the problem of bound inhomogeneities to be 
addressed more directly without the complication of gravitational wave inhomogeneities. 

The choice of spherically symmetric inhomogeneity does not permit exact solutions of Einstein's equations 
except where fluids have vanishing pressure. However, a clear physical picture of the behaviour of spherically 
symmetric inhomogeneities with non-zero pressure can be obtained by the device of matching together 
homogeneous solutions of different curvature and density. The resulting patched solution describes the 
evolution of spherical overdensities (or underdensities) in a smooth background universe. In the limit that 
the inhomogeneities become small, they will evolve in accord with the results of small perturbation theory. 
When the inhomogeneities are not small, we obtain a description of non-linear processes like void formation, 
the condensation of gravitationally bound lumps, or the creation of primordial black holes. This technique 
was first introduced into general relativity to study the evolution of inhomogeneities by Lemaitre [|ll]] and 
was subsequently applied to the study of protogalaxies by Harrison |12[| . Here, we shall apply it to the 
equations of string cosmology to further our understanding of the evolution of the pre- and post- big-bang 
phases in the presence of spherical inhomogeneities. 

The string cosmological models considered here are derived from the bosonic sector of heterotic string 
theory reduced to (3+1) dimensions of spacetime. They are assumed to have vanishing cosmological constant 
and vanishing Maxwell field. Their field content consists of an antisymmetric tensor field, a dilaton, and the 
space-time metric tensor. However, in the low-energy limit only the antisymmetric tensor field strength is 
important in the equations of motion. 

In 10-dimensional superstring theory, gravitational anomalies occur which signal the breakdown of 
energy-momentum conservation. However, by redefining the antisymmetric tensor field strength, these 
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anomalies can be cancelled. The redefined antisymmetric tensor field strength, H, is given by [13] 



H = dB + uj l , 

where B is the antisymmetric tensor field and ujl is the Lorentz-Chern-Simons (LCS) form involving the 
Lorentz spin connection. The antisymmetric tensor field strength is a 3-form. In four dimensions it is dual 
(in the sense of differential forms) to a one-form which can be shown to be the gradient of a scalar field, the 



axion. FRW spaces are maximally symmetric and so a theorem proved in [14] implies that all contributions 
from the LCS terms vanish. Hence, the energy-momentum tensor for the axion and the dilaton in the 
Einstein frame consists of two coupled stiff perfect fluids. 

The low-energy limit of string theory provides a new picture for the evolution of the early universe. 
Once a stage of low coupling and small curvature is reached, the universe enters the 'pre-big-bang' era [15]. 



During this stage the universe undergoes super inflation (accelerated expansion) in the string frame driven 
by the kinetic energy of the dilaton. By contrast, in the Einstein frame this corresponds to an accelerated 
contraction. The pre-big-bang era ends when the string coupling becomes strong enough for the low-energy 
limit to be no longer valid. However, exactly how this "graceful exit" can be effected is not yet fully 
understood [pUf !. In the transition era, complicated non-perturbative effects will become important and it is 
not clear if a curvature singularity can always be prevented by higher-order contributions. Eventually, the 
universe must enter the (classical) post-big-bang era. Therefore, the pre-big-bang phase can be understood as 
a way to provide initial conditions for the classical post-big-bang era. An interesting aspect of this scenario 
is provided by the duality symmetries present in string theory. In spatially homogeneous cosmological 
models, scale-factor duality relates solutions for the pre-big-bang phase ("+ branch") to those for the post- 
big bang phase ("- branch"). Unless there is a self-dual solution there is the problem of how to relate 
these two branches. This may require an explicitly quantum cosmological transition (T^j. In the presence 
of inhomogeneities, especially those which allow some parts of the Universe to expand whilst other parts 
collapse, the impact of duality invariance may prove more unusual and motivates further detailed study of 
realistic inhomogeneous string cosmologies. 

Here, we extend our understanding of inhomogeneous string cosmologies by investigating the simple 
model of non-linear spherically symmetric inhomogeneities outlined above, in which a spherical curvature 
perturbation is self-modelled by an FRW universe of non-zero curvature in a flat background FRW universe. 
In section 2 we describe the self-modelling of spherical inhomogeneities in general relativistic universes 
containing fluids with pressure equal to density. In section 3 the connection with string theory is displayed. 
In sections 4 and 5 the post and pre-big-bang solutions are given and, finally, the results are discussed in 
section 6. 

2 Spherical Inhomogeneities in General Relativity 

Spherically symmetric density inhomogeneities can be modelled by matching a section of a closed (or open) 
FRW universe to a background universe described by a flat FRW universe in such a way that the metric 
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and its derivatives are continuous at the boundary. First, consider this matching in general relativity for the 
simple case of a universe containing a perfect fluid, with pressure p and energy density p, which are related 
by a stiff equation of state, p = p. 

The flat background FRW universe has an expansion scale factor R(t) and its dynamics are determined 
by the Friedmann equation 
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dt J 



SV = ^i? (i) 



where t is proper time. Note that Einstein's equations are invariant under time-reversal, i.e. R(t) = R(—t). 
Energy-momentum conservation for the p = p fluid implies that 

pouBT 6 . (2) 



and so R(t) oc £3. 

The dynamics of the perturbed region of non-zero curvature are described by another Friedmann equa- 
tion, with a scale-factor S(r) 

dS\ 2 8-irG, , 



dr 



(p + Sp)S 2 - k, (3) 



where r is the proper time within the density perturbation, dp, and the constant, k, measures its spatial 
curvature. 

The proper time, r, inside the perturbation and that in the background universe, t, can be related using 
the equation of relativistic hydrostatic equilibrium |l2| ]l8| 



dr p + p 

where <&( 9rav ) is the (Newtonian) gravitational potential, and r is the radial distance, and so 



dr = exp[& 9rav) ]dt. (5) 

The equation for hydrostatic equilibrium for a perfect fluid is derived under the assumptions that the 
configuration is static (that is, only spatial derivatives are non-vanishing) and that the gravitational field 
is weak, so that the Newtonian gravitational potential &(9 rav ) completely determines the metric. In this 
case equation (Q) follows directly from the conservation of energy-momentum for a perfect fluid. In this 
particular matching problem of two Friedmann universes staticity means that there should be no radial flow 
of matter in or out of the perturbation. 
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Equations @ and (|5|) imply 

dr /S^ 3 



where we have introduced the constant 5, the density contrast parameter, defined by 

5=*A I*I<1- 
p 

For the case of an over density we have 5 > 0. 

Assume that at some initial time to the perturbation appears and the following matching conditions hold 
between the metric scale factors inside and outside the perturbation, jT^ , 

*- ©„-(f)„- 



Then, we have 



dRy G 
~dt) ~R*' 



(8) 



^(P-ES% (9) 



- 



dr J 

where we have defined three new constants by 

6 K 
Furthermore, this implies that S is given in terms of the background scale-factor, R, by 

where we have defined a further constant by 

F ^ Hg(l + 6b)" 

Equation (^implies that the maximum of the scale- factor of the fluctuation is S+ = F~ 1 , and hence 
the integral of ( |10|) is given by 

s = s + ^Cf/f) (12) 
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where R* is the value of R at S + . The fluctuation thus begins expanding with the background but is slowed 
with respect to it because of its overdensity. Eventually, its expansion is halted by its self-gravity and it 
begins to collapse whilst the background continues to expand. 

We see that the spherical fluctuation has vanishing scale-factor S, i.e. zero radius, for R = and also as 
R — » oo, where S oc R . This means that the perturbation does not collapse to a black hole in finite proper 
time, since R oc is. Physically, this situation arises because the p = p fluid possesses a Jeans length equal 
to the horizon size and the fluid pressure is able to resist gravitational collapse as soon as the fluctuation 
enters the horizon. This situation is distinctive compared to the evolution of overdensities when p < p. 
In the p < p case significant overdensities collapse to form black holes if they turn-around when they are 
intermediate in scale between the particle horizon and the Jeans scale |20|]. In practice, even in the p = p 
case some black-hole formation might be possible in finite time because of small fluctuations in the sound 
speed and the horizon size but, realistically, we would expect shock formation to play a role in damping 
non-linear inhomogeneities in the p < p cases (see also the discussions of the scaling properties of black hole 



formation found by Choptuik and others [21] in this connection). 

Our description of the evolution of the spherical overdensity is completed by deriving the ratio of the 
density in the perturbation to that in the background universe. This is given by 

p (pert) {R ' 

-jfcS) =(1 + 5o) [s 

Using equation (|T2|), this can also be expressed as 



7^ = (1 + <5o) VMs, 



i — 

R* 



1 3 

(13) 



If we take the limit of small time (or small R) then we recover the usual description of the growth of 
small perturbations in time in an appropriate gauge. 

Finally, it is interesting to note a simple duality scaling property that appears in the above analysis. 
To see it in context, we can generalise the analysis given above to the case of inhomogeneities in a general 
perfect fluid model with equation of state p = (7 — l)p. Equation (|i~0|) then generalises to 



AR) \R> a+W^H-^J 1 («) 

We see that only in the case of a stiff perfect fluid, that is 7 = 2, does ( |i~4l) admit both a scaling symmetry 
(S — > aS, R — > aR, a constant) and the duality invariance 

R — > R" 1 . (15) 

In Figure 1 we show the evolution of R, R , and S with respect to the background proper time, t. The 
background expands as a flat FRW universe (R oc t 1 ^ 3 ), whilst the overdensity expands less rapidly, reaches 
an expansion maximum, but then collapses more slowly, tending to zero size at an infinite future time. 
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3 Connection with String Cosmology 



Since FRW universes are spherically symmetric about each point, the LCS contributions to the antisymmetric 
tensor field strength vanish, and the dilaton and axion fields behave as two coupled stiff (7 = 2) perfect 
fluids in the Einstein frame. The Einstein frame is related to the string frame by a conformal transformation 
of the metric of the form, 

^Einstein ^String 

9^ v = e g^ v y , (lb) 

where $ is the dilaton. The advantage of the Einstein frame is that it is the frame in which the action takes 
the Einstein-Hilbert form and so, if there are no LCS contributions, the problem is that of general relativity 
with an energy-momentum tensor consisting of two coupled stiff perfect fluids. However, in order to provide 
a physical interpretation of the dynamics, we should go back to the string frame. In string theory the 
equation of motion for a classical test string implies that its world sheet is a minimum surface with respect 
to the string metric; that is, the (two-dimensional) analog of a geodesic for a point particle. Therefore, it 
appears that strings "see" the string metric rather than the Einstein metric j22|. 

The low-energy effective action in the Einstein frame yields to the following set of equations (k 2 = 
8ttG,c = 1) 

R^-\g^R = ^(Wt^+^t^) (17) 

V^e-^H^) = (18) 
6' 

where 



□$ + i e - 2 *i^ A #^A = (19) 



^T,u = \{$,^, u -\g» u m 2 ) (20) 
{H %u = ^e- 2 *(3H, XK H u ^-± 9flu H a(37 H^) (21) 

In four dimensions, if we take the space-time dual of the antisymmetric tensor field strength H^ u \, the 
Bianchi identity dH = shows that the dynamical content of the antisymmetric tensor field strength is 
determined by a (pseudo-) scalar field, namely the axion field b. So, we may write 

where e^ AK is the totally antisymmetric Levi-Civita symbol. Furthermore, in the pure dilaton model we 
will discuss, the perfect fluid is characterised by equal density and pressure, p = p = |$ 2 , and the 4-velocity 
obeys u a u a = — 1. 
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Solutions for flat and closed FRW models have been discussed in 0. Using conformal time rj the metric 
for a FRW model with a scale factor a{rj) can be written as 

ds 2 = a^^i-dri 2 + + r 2 dQ 2 ). (22) 

1 — kr z 

Naturally, the general-relativistic scenario discussed in the previous section translates directly into a 
description of the post-big-bang era. However, since the pre-big-bang model makes use of the invariance 
of solutions under time reversal, it is misleading to make reference to S+ and Rq = So since these are two 
"time-ordered" length scales. If initial conditions are specified at Rq = So and developed into a state at S + 
then the time reversal causes 5+ to lie before Ro = So; that is, conditions at Ro can no longer be treated 
as initial conditions. Therefore, in the discussion of pre-big-bang solutions only the initial scale -Ro = 5o 
should be used in the solutions for the scale factors. 



4 Post-Big-Bang Solutions 

For simplicity, only overdense perturbations {5p > 0, k > 0) will be discussed in this section. The symbol r 
denotes conformal time in the flat background universe while 

r (in) ig 

the conformal time inside the spherical 
perturbation. Furthermore, the index "*" refers to the epoch when the scale factor of the perturbation S 
has its maximum value, 5+. Proper time runs from — > +oo. Solutions for a pure dilaton model are given 
in §• 

(i) Outside the perturbation, modelled as a flat FRW universe: 
The scale factor is given by 

R 2 (t) = Rl- (23) 

T* 

and the dilaton evolves as 



e*= L) (24) 



(ii) Inside the perturbation, modelled as a closed FRW universe 

Introducing a function rj of the conformal time r( m ) inside the perturbation defined by 



, ? (r( in )) =tan(r( in >), 

the scale factor is given by 



S 2 (r 1 ) = 2S 2 + -^— 2 (25) 
1 + rj z 
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and the dilaton behaves as 

e* = r] ±VE , (26) 

since 77* = 1. 

In order to relate the conformal times inside and outside the fluctuation, one can make use of equation 



(12), which gives 



- = V . (27) 

Furthermore, the evolution of the background scale-factor in proper time t, in the Einstein frame, is found 
after integrating dt = Rdr, to be 

1 

t \ 3 



R = R* y-j (28) 

As discussed above, the frame appropriate for physical interpretation is the string frame. In this frame, 
the background scale factor is given by 

(S)R = e */2 R = R„ (L\ . (29) 

Expressed in proper time in the string frame, i.e. after integrating d( S 't = ( S 'Rdr, ( s 'R(( s >t) is given by 

iS,R = R -{m) ■ < 30 » 

Pre- and post-big-bang, respectively, refer to the background universe. In the post-big-bang epoch the 
upper sign in ( |30| ) is chosen. The matching conditions (cf (|7|)) provide a well-behaved metric in the Einstein 
frame. Since physics should not be frame dependent, it is expected that similar matching conditions should 
hold in the string frame. As one can show, this implies that the sign chosen in ( |26| ) for the evolution of the 
dilaton inside the perturbation should be the same as that in the background universe; hence 



exp (V 6 ^)) = exp (Vf ert ) 



(31) 



The energy densities in the Einstein frame and the string frame can be related by going back to the 
definition of the energy-momentum tensor as a variational derivative of the Lagrangian C [23]. In the 
Einstein frame, we have 
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-too 



HO" 



5C 
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Using the conformal transformation given by equation (jig), we find that the energy density in the string 
frame is given by 



(S) -2$ 
K * p = e p. 



(32) 



Hence, using fl3T|), we confirm that the ratio of the energy density in the perturbation to that in the 
background is the same in the Einstein and the string frames: 



(S)p(pert) p{pert) 
(S) p{back) p(back) ' 



(33) 



5 Pre-Big-Bang Solutions 

In order to find solutions depending on one length scale we begin again with equation (|l~0|). Furthermore, 
since the analysis does not depend on a maximal length scale, the solutions for underdense regions (modelled 
by matching to an open FRW universe, with k < 0), can be treated as well. 

5.1 Underdense Regions 

In this case the density parameter 5 is negative and hence the constant F, defined in (11), is also negative. 
Equation (|lO|) now becomes 



dS__S_ 
dR~ R 



1+\F\S 4 V . 



(34) 



This integrates to give 



S_ 

s 



^[(l + *or5-l]*(i) 



\So\ 



l+<5o 



(l + <5 )"*-l 



2 ^ 4 
Ro 



(35) 



Hence, the ratio of the density in the fluctuation to that in the background (in the Einstein frame) is 
given by 



P 



(pert) 



1+^0 



J back) 



(l + 6o)-? -1 



l + S 



ri + 6o)~3 -i 



( R_ x 1 



(36) 
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5.2 Overdense Regions 

In this case both the density parameter, 5, and the constant F are positive. An integration of equation ( |l0|) 
yields 



V2 



1 - (1 + 6 



(ro) 



So 



1+So 



1 - (1 + S 



(37) 



and the ratio of the densities in the Einstein frame is given by 



9 (pert) 
Jback) 



l + So 



1 - (1 + So 



So 



1 + 5, 



+ 



1 - (1 + S 



R - 1 



(38) 



In order to connect the conformal times inside and outside the perturbation, and hence the solutions for 
the pure dilaton model, we will use the original notation of ||. 
(i) In the flat background the scale factor is given by 0] 

R 2 = ^=t (39) 
^3 

with K some constant and r the conformal time in the background. The evolution of the dilaton is given 
by 



..= £j ,40) 

where r/ is an integration constant. For pre-big-bang solutions the - sign is chosen. In terms of proper time 
in the Einstein frame these solutions read 



±V3 



2 

t \ ^ 



and 

R 2 = ( 

2V3 



R 2 =[~t)\ (42) 
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Recalling that time is negative in the pre-big-bang era, one sees that in the Einstein frame the background 
universe is contracting and the dilaton is described by a growing function. The scale factor in the string 
frame is given by 

^R 2 = e*R 2 = Q) § tf (^L) § tl^\ (43) 

and hence the universe is expanding in the string frame. 

(ii) Inside the density fluctuation the scale factor is given by 



,2 K 



with K some constant, and 

V 



= — 1 (M) 



tan(rj n ), k = +1 overdensity 
tanh(rj n ), k = — 1 underdensity 



where Tj n is the conformal time inside the fluctuation. Using the matching conditions, ([?]), tq and rjo can be 
related at this epoch by 



The solution for the dilaton is found to be 

X -y/3 




(46) 



with rjf an integration constant. 

Since the dilaton solutions are involved in the transformation of the energy densities from the Einstein 
frame to the string frame, it is enough to relate r and r\ (and not explicitly the conformal times inside and 
outside the fluctuation) and this can be done using equations (|35|) and (|37[). First, note that 



and 



S\ 2 _ yjl + kyj) 

So) vo(l + W) ■ { ' 
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Observing that 5q is positive or negative, according to the sign of k, equations ( |35[ ) and (37) can be 
reduced to a single expression, namely 



b \ w - T (L) (49) 



Ro 



7 + «M 7ET 

with 



a=l-(l + <y )~' , 7= 



i 



1 + So 
using ei 

expressed in terms of the conformal time r in the background universe, 



The function T(^) is defined in accord with equation (|47D . So, using equation (f48|), the function 77 can be 



1 - Jl - 4kE2T 2 
V = k V ^ , (50) 



where 

■p — 

an = 



1 + kril 

which can be expressed in terms of tq using ([45]) . 

Therefore, in this more general case, the relationship between the dilaton solutions inside and outside 
the fluctuation is more complicated than that found earlier for the post-big-bang era (cf equation (|3l"|)). 

The perturbation originates at some epoch tq with a scale factor corresponding to a value Ro in the 
background universe. The big-bang occurs at r = 0, but one should not follow the evolution all the way 
to r = since at some point the string coupling becomes too strong for the low-energy action to remain 
a valid approximation to the full theory. In the pre-big-bang phase the proper time is restricted to t < 0. 
Suppose that a perturbation originates at to and the period of interest is between to and some time t s , (with 
I t s \<\ to |), when the string coupling becomes strong, so we are interested in the regime t < tq, R < Rq. 
Assuming R <C Ro, the scale factor of the perturbation S and the function of its conformal time 77(7") can 
be expanded as follows, 

S /2a\v R 



so - v 7 ; * (51) 
^ * ^(i + ^r 1 -- (52) 

Vo 7 T o 



Thus, ^ ~ -j^, ^ ~ ^. However, for <5o > 0, the proportionality factor in ( |5l| ) is bigger then unity, and 
so the scale factor of the perturbation exceeds that of the background universe (see figure 2). 

This might be interpreted as interchanging the role of the flat and the closed FRW models (ie there 
are inhomogeneities corresponding to flat space-time sections in a closed FRW background universe ||24|). 
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Figures 3 and 4 show the behaviour of the two scale factors, R and S, with time, t, for different choices of 
the initial density parameter 5q. In Figure 3, the scale factor of the overdense region is larger than that of 
the background universe within the period of interest between to U P to some time, t s . In Figure 4, where 
there is an initial underdensity at to, the scale factor of the perturbation stays below that of the background 
universe. 

It is clear that the evolution of the two scale factors is very similar, from the origin of the perturbation 
up to the big bang. Thus, one cannot actually say that the perturbation evolves independently from the 
background. Qualitatively, the same behaviour is found in the string frame. The evolution of the dilaton 
depends on the respective conformal times (or functions of them), namely r and r\. However, in the regime 
between the origin of the perturbation and the big bang these are directly proportional to each other as can 
be seen from (|5^). Hence, the evolution of the scale factors in the string frame is very similar as well. 

In accord with the behaviour in the general relativistic (hence post-big bang) regime, no primordial 
black holes should be formed in the pre-big-bang era in this model since the perturbation does not evolve 
independently and so does not form trapped surfaces during the collapse of the scale factor of the background 
universe in the Einstein frame. 

From this discussion it can be concluded that the introduction of a spherical perturbation modelled as a 
closed or open FRW model in a flat FRW model does not destroy the global isotropy of the flat background 
universe. Hence the pre-big-bang flat FRW solution is robust with respect to this special type of curvature 
perturbation. 

In heterotic string theory, any solutions containing only a dilaton can be transformed into solutions 
containing a dilaton and axion without changing the background metric. S-duality reflects the fact that 
heterotic string theory is invariant under SL(2, IR) transformations. Define a complex scalar A = b + ie - * 
then g§ 

\ V = + f gn U -> gnu a, f3, 7, S e IR, a5 - f3<y = 1 

is a solution to the low-energy equations of motion (equations (|T^)-(|l9|)). 
If we start with a pure dilaton solution 

A = te 

and choose a = /? = 7 = — 5 = l/y/2, then A' = b new + iexp[— & new ] is given by 

b new = tanh$ (53) 
e * n ™ = cosh*. (54) 



Hence, the pure FRW dilaton solution has 
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where r is a different function of the conformal time r\ for each of the FRW models which can be read of 
from the solutions given above, and is related to a FRW dilaton-axion solution with 




which is given in Q. 

However, near the big-bang (that is, for small ^), the dilaton solution goes over into that obtained 
for the pure dilaton case and so the conclusions from the pure dilaton case are not likely to be changed 
fundamentally. 

6 Conclusions 

A very simple inhomogeneous string cosmological model has been investigated in the pre- and post-big- 
bang eras. In the post-big-bang era the usual general relativistic behaviours of a growing overdensity or 
underdensity are found. This leads to an increasingly inhomogeneous universe if the initial inhomogeneities 
are significant, or they are not inflated away. However, in the pre-big-bang phase the global isotropy of the 
background FRW model is unaffected by the introduction of a spherically symmetric curvature perturbation. 
This rather different behaviour of the model in the pre- and post-big-bang stages is a manifestation of 
the fact that these two regimes pick out different parts of the general relativistic solution. In general 
relativity, considering only positive (proper) times, the solution is naturally divided by the time to when 
the perturbation originates. The physical solution is then given for times later than to, whereas the part 
between the big bang and to is discarded. However, in the pre-big-bang stage one is considering a time- 
reflection of the general relativistic solution. Hence, the part of the solution that was discarded in the 
general relativistic case (post-big-bang) becomes the physical solution in the pre-big-bang era. Thus, it is 
the physical solution of general relativity which is considered unphysical in the pre-big-bang regime. Hence, 
unlike in the general relativistic solution, a curvature perturbation like that considered here does not lead 
to an inhomogeneization of the universe in the pre-big-bang epoch. 

The authors of references [|| and || discussed inhomogeneous pre-big-bang models using standard ap- 
proximation techniques of general relativistic cosmology. They found that the approximation of neglecting 
spatial gradients becomes better as the big-bang singularity is approached. This conclusion is confirmed by 
the simple quasi-homogeneous model discussed here. The evolution of a curvature perturbation modelled 
by a section of a Friedmann universe is very similar to that of the background Friedmann universe and 
indeed becomes more similar to it with time. Hence isotropy is conserved and spatial gradients do not be- 
come important. We note also that part of the general solution for general relativistic cosmological models 
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containing a p = p perfect fluid is a perturbation of a Kasnerlike solution which (unlike for case of perfect 
fluids with equation of state p < p) contains the isotropic expansion as a particular case. Because of this, 
the general behaviour of the equations of string cosmology at early times is significantly simpler than that 
of general relativity in vacuum. 
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Figure 1: Evolution of the background scale factor R oc ta , the transformed one R 1 oc t a and the scale 
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factor of the perturbation S(t) oc ts (1 + £s)~2 in terms of proper background time t. 




Figure 2: The proportionality factor ^2^ J as a function of 5q- 
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Figure 3: Evolution of the scale factors of the perturbation S and the background universe R as a function 
of proper time t in the Einstein frame ; 5q = 0.5, Rq = 1, to = —1. 




Figure 4: Evolution of the scale factors of the perturbation S and the background universe R as a function 
of proper time t in the Einstein frame; 5q = —0.5, Rq = 1, to = —1. 
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